Effect on Cavity Optomechanics of the Interaction Between a Cavity Field and a ID 

Interacting Bosonic Gas 



Qing Xing-Hua Ru\ W. M. Uu\ X. C. Xie^'"^, and An-Chun Ji^'i 

^Beijing National Laboratory for Condensed Matter Physics, 
Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 
^Department of Physics, Capital Normal University, Beijing IOOO48, China 
^International Center for Quantum Materials, Peking University, Beijing 100871, China and 
^Department of Physics, Oklahoma State University, Stillwater, Oklahoma 74078, USA 

(Dated: January 15, 2013) 

We investigate the optomechanical couphng between ID interacting bosons and the electromag- 
netic field in a high- finesse optical cavity. We show that by tuning the interatomic interactions, one 
can realize the effective optomechanics with the mechanical resonators ranging from the side-mode 
excitations of a Bose-Einstein condensate (BEC) to particle-hole excitations of Tonks- Girardeau 
(TG) gas. We propose that, this unique feature can be formulated to detect the BEC-TG gas 
crossover and measure the sine- Gordon transition continuously and nondestructively, which are 
achievable immediately in current experiments. 
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The experimental achievements in manipulating the 
strong coupling between ultracold atoms and the electro- 
magnetic field in an optical cavity have triggered many 
new exciting advances to cavity quantum electrodynam- 
ics (QED) [iHsl. One of the remarkable achievements is 
to implement cavity optomechanics with cold atoms [gI, 
or a BEC [sl] , which is of great importance both for tech- 
nical applications ranging from optical communication to 
quantum computation [9], and conceptional exploration 
of the classic-quantum boundaries [10| . 

In this Letter, we investigate the optomechanical cou- 
pling between a ID interacting bosonic gas and a cavity 
field. So far the recent works have neglected the inter- 
atomic interactions or considered merely the weakly in- 
teracting region, where the mean-field Bogoliubov theory 
is valid [8|, luj. In this case, the ID bosonic gas forms 
a BEC (or quasi-condensate) [l2[, and the side- mode 
excitations of the condensate play the role of mechan- 
ical resonator with the bare frequency uj^ = Ahk'^/M 
0, here k = 27r/Ac is the wave- vector of the cavity 
mode. However, when the interatomic interactions are 
added into the system, the situation changes dramati- 
cally. The strong interatomic interactions would trans- 
form the ground state of the condensate to a Luttinger 
liquid (LL); and remarkably in the strongly interacting 
limit, the ID bosons-known as a TG gas jl3l4l6| -exhibit 
completely different behavior like ideal fermions. It is 
therefore important to explore the interatomic interac- 
tion effects on the cavity optomechanics, where the quan- 
tum fluctuations of ID bosons are very strong. 

In this work, we first employ the quantum hydrody- 
namical approach to derive an effective model of the 
cavity QED with ID interacting bosons. We show that 
the effective optomechanics can be realized in the inter- 
mediately and strongly interacting regions. The corre- 
sponding optomechanical coupling is determined by the 



low-energy excitations or the interatomic interactions of 
the bosons. Therefore, by probing the cavity oscillations 
or the noise spectra versus the interatomic interactions, 
one can determine the quantum phases of the ID in- 
teracting gas and detect the BEC to TG gas crossover, 
a fascinating phenomenon of the system [17|. Further- 
more, we propose that one could also measure the sine- 
Gordon transition, which has stimulated considerable in- 
terest 18|, I19L conveniently with the nondemolition mea- 
surements [20|, |21|. 

The system under investigation is schematically de- 
picted in Fig. [IJa), where N ultracold bosonic atoms 
of mass M with resonant frequency uJa are confined in 
a ID trap inside an optical cavity with length L. The 
cavity mode of frequency uJc is driven by a pump laser of 
frequency ujp at rate 7^, and k is the decay rate of the cav- 
ity field. Following Ref. jlH, we adiabatically eliminate 
the internal excited state of the atoms, as justified by the 
large detuning between the atomic resonance and pump 
frequency. Then, by using the dipole and rotating-wave 
approximations, one arrives at the following Hamiltonian 
of the atomic part 
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Here, ^(x) is the bosonic field operator and V{x) = 
hUo cos^{kx)c^ c is the dynamical periodic potential, with 
c the annihilation operator of a cavity photon, and the 
potential depth Uq = g^/iuOp — oOa)- The interatomic in- 
teractions are given by contact pseudo-potentials U{x — 
x') = giJix - x'), where gid = (i-Ca.TX)M/i 
the effective ID coupling strength with ag the three- 
dimensional scattering length, C = 1.0325, and l± = 
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FIG. 1: (color online). Experimental set-up and schematic 
density distribution, (a) N bosonic atoms with resonant fre- 
quency oja are confined in an effectively ID trap inside an 
optical cavity of length L. The cavity mode is driven by 
a pump laser of frequency ojp, and k is the decay rate, (b) 
Ground-state atomic density distribution po{x) in the absence 
of a cavity mode, (c) Distribution of density fluctuation pq{x) 
with wave- vector q = ±2/c, which is scattered by the periodic 
potential of the cavity mode. 



^JhjMujx. the transverse osciflator length. 

We start by considering the general situation with ar- 
bitrary interatomic interactions and derive an effective 
model of the the system by using quantum hydrodynami- 
cal approach [23] , which is a well-defined low-energy the- 
ory. We shall work in the low photon numbers limit, 
where the dynamical periodic potential V{x) is feeble. 
By introducing two new fields (\){x) and B{x\ which de- 
scribe the collective fluctuations of the phase and den- 
sity respectively and satisfy the commutation relation 
[^(x), ^dx'^{x'y\ = i^[x — x'^^ we can express the bosonic 
field operator ^(x) as 
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and the corresponding density operator 
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In the following, we shall first consider the long- 
wavelength approximation, i.e. Ac ^ 1/po, where we 
can only keep the m = term in Eqs. ^ and ([3]). In 
this limit, the system can be expressed by the hydrody- 
namical description, with the weak dynamical periodic 
potential coupled to the slow part of the density oper- 
ator p{x). Then, the corresponding low-energy effective 
Hamiltonian of the atomic part reads 



1 



H'a = 1^ dx{'^[K{dJ{x)f^-{dJ{x)-npo)'] 



V{x) 



d,0{x)}. 



(4) 



Here, K is the dimensionless parameter and Vg the sound 
velocity. They both depend on a single dimensionless in- 
teracting parameter 7 = Mgid/h'^po, with VsK = vp = 



hnpo/M fixed by Galilean invariance. We note that, the 
Hamiltonian (j4]) describes a LL coupled to a weak pe- 
riodic potential, which is dynamically dependent on the 
atomic state and determined self-consistently. 

We further transform the Hamiltonian (|4]) to momen- 
tum representation and then implement the standard 
bosonization procedure by introducing the bosonic cre- 
ation operator 6j = '\J -^^Pq- We can arrive at the fol- 
lowing effective Hamiltonian of the coupled system 
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Here the first term describes the long- wavelength density- 
fluctuations of the ID interacting gas with Uq = Vs\q\ for 
\q\ <C Pq^. The second term is the coupling between the 
corresponding density-fluctuations and cavity field with 

g = ^ In the sums, we assume that only the q = 

±2k modes are coupled to the cavity, which is justified 
by the low photon numbers limit. A = Uc — ujp -\- UqN/2 
is the effective cavity detuning. 

The above effective Hamiltonian (|5]) actually describes 
the optomechanical coupling between a mechanical oscil- 
lator (with the frequency uom = ^±2k = 2kvs) and the 
radiation pressure force of a cavity field. To see this, 
we introduce the quadratures of the bosonic excitations 
= Z]g=±2/c(^J + ^g)/V^7 and derive the following 
Heisenberg-Langevin equations 
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with the resonance frequency A eff = A + \/2gXM and 
V^Cin the noise term. Here, the low-energy long- 
wavelength density fluctuation (phonon) plays the role 
of a mechanical resonator. 

Now, we investigate the optomechanical coupling gov- 
erned by the set of coupled Eqs. (|6])-([7]). One of the 
characteristic phenomenon of cavity optomechanics is the 
bistable behavior of the stationary solutions, which we 
derive as Xm = —2^/2g\c\^ / uj m -> and |cp = rf /[ti^ + 
(A - Ag'^uj^\c\^Y]. In Fig. [2^, we give the bistability 
for three typical oscillator frequency cjm, which corre- 
spond to different quantum phases of the system (see be- 
low). The linear stability analysis shows that the middle 
branch (dashed line) is unstable, while the up and low 
branches are stable. Here, the typical experimental pa- 
rameters are used: L ~ 100 /im, Ac = 780 nm, N c:^ 5000 
^^Rb atoms, with 1/po — 20 nm<C Ac satisfing the wave- 
length approximation; = 27r x 1 MHz, /7o = 27r x 20 
kHz with the Rabi frequency go = 27r x 10.9 MHz and 
the pump- atom detuning cjp — 0;^ = 27r x 32 GHz. 
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FIG. 2: (color online). Steady- state and dynamical 
behavior of the optomechanical coupling. (a) Mean 
cavity photon number versus A for ujm = 27v x 
[0.02 (BEG), 0.12 (LL),0.28 (TG gas)] MHz. Here r/cr is the 
bist ability threshold, (b) Gavity photon number |q;|^, coupled 
oscillator Xm (dashed line), and free oscillator XM{g = 0) for 
A = 0. (c) Noise spectrum Sxc versus cj and A. (d) Gen- 
ter frequency ojs of versus 7. Inset shows the results for 
small 7 region. 

To discuss the dynamics of the optomechanics, we have 
taken into account the lossy and driven cavity, where 
quantum jumps in the cavity photon number can lead to 
a strong entanglement between the cavity photon number 
and bosonic wavefunction. This creates a displacement 
noise spectrum of the mechanical oscillator Sxm{^) — 
2/^(4^ccc;m)^[/^^ + (A + a;)^]/|(i(cc;)P, and the correspond- 
ing measurable noise spectrum of the cavity field quadra- 
ture Xc = (ct + c)/V2 0, which gives SxA^) = 
[(2^cA)25xM(^)+2A.[/.2+(A+c^)2]]/|^i»|2 (see Fig. Efc). 
Here, d{uj) = (cj^ - cj|^)[(/^ - i^Y + A^] -f 2cjmA(2^c)^, 
and A = A - (2gcf jijju- 

We then explore the interatomic interactions depen- 
dence of the effective optomechanics. First we note that, 
whether the system is in weakly or strongly interacting 
regions, one can realize the effective mechanical oscil- 
lator in the whole regions. This implies that the ID 
bosonic gas is in fact in a universal class, which can be 
well described by the low-energy hydrodynamical the- 
ory. However, the mechanical oscillator frequency cjm is 
fully dependent on interatomic interaction parameter 7 
with um = 2/c^;5. Here, the sound velocity is given by 
Vs = V(po/M) d^E/dN^, and E = Y.i h^kf/2M is the 
ground-state energy of ID bosonic gas with ki determined 
by the following Bethe-Ansatz equations [1^ 

kiL = 2^Ii - V tan-i , (g) 

where li G {-(A^ - l)/2, • • • , (TV - l)/2} are the set of 
integers. We solved Eq. (|8]), and in Fig. [3] we show the 



numerical results of ujm versus 7. Experimentally, when 
photons enter the cavity, the light filed and the bosonic 
excitations are coupled nonlinearly, the eigenfrequency of 
Eqs. (6)- (7) will modify uw Nevertheless, the optome- 
chanical coupling is fully determined by the frequency 
cjm of bosonic excitations and thereby dependent on 7 
of the ID gas. Accordingly, by measuring on the cavity 
field oscillations or detecting the noise spectrum, one can 
determine the continuous BEC-TG gas crossover, which 
is an intriguing result of the system. 

In the weakly interacting region, we can use 
the B ogoliubov ap proximation to derive ujm = 
2kvF\/^ — 7^/^/(27r)/7r for 7 < 10 (the dashed line of 
Fig. [3]). However, we note that uom vanishes as 7 ^ 0, 
which contradicts with recent experiments (sf. In fact, 
as we know because the bosons begin to (quasi) con- 
dense for 7 <C 1, the dominant contribution of the den- 
sity fluctuation p±2k = ^±2fe^^=o + Eg/o ^±2^+^^^? is 
the quasi-particle excitations from macroscopic occupied 
q = ground state. The energy of a quasiparticle is 
then determined by the Bogoliubov excitation spectrum 
um = \/e±2/c(e±2/c + 2^i(iPo)/^ (the dotted line of Fig. 
[3]). In this case, the mechanical oscillators are the side- 
mode excitations of a BEC; and in the limit of 7 = we 
get the bare oscillator frequency uj% = Ahk'^ /M (inset of 
Fig 13]). Accordingly, the coupling between the oscillator 
and a cavity field should be replaced by ^ = N /2 

in Eqs. (|6j)-([7]), which is enhanced by the condensation. 
We give the stationary bistability in the inset of Fig. [2^, 
and then numerically integrate the coupled equations for 
A = by switching on t^//^ = 4 at t = 0. We find that 
the cavity field oscillates regularly (not shown), but the 
frequency of the cavity field has a large shift of uom , which 
agrees well with Ref. (sf . 

When 7 is further increased above unity, the collec- 
tive density excitations become dominant and the ID gas 
crosses to a LL |17]. Numerical integration of Eqs. ([6])- 
([7j) shows that both \a\^ and Xm exhibit well-defined 
oscillations, and the cavity field is excited resonantly at 
Xm = (Fig. Wp)' Yet we note that, different from the 
BEC phase where g is collectively enhanced, the optome- 
chanical coupling becomes small. In this case, the linear 
stability analysis shows that the eigenfrequency of the 
set of coupled equations is nearly the same with the free 
oscillator cjm, see Fig. [JJd for example. This is a char- 
acteristic phenomenon of the region. Experimentally, by 
increasing the interatomic interactions, if the oscillation 
frequency of the cavity field follows the the solid line of 
Fig. [31 the ID gas should be in the LL phase. We also 
calculate the noise spectrum 6'xe (Fig. [2t), where the 
center frequency uos of the spectrum has a shift of uom- 
In Fig. [2]i, we give uos for A = in the whole interacting 
regions. We find that uos increases with 7 and has a ten- 
dency to saturate above 7 50, which can be inspected 
in experiments. 

While for 7 ^ 1, the strong interactions would 
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FIG. 3: (color online). Mechanical oscillator frequency ujm 
versus 7. The solid line is the numerical Bethe-Ansatz re- 
sult. The dashed and the dash-dotted lines are the asymp- 
totic results in the weakly and strongly interacting regions 
respectively. In the weakly interacting region, we also give 
the Bogoliubov excitation spectrum (the dotted line). 



prevent the bosons from occupying the same position. 
Especially when 7 = 00, the symmetric many-body 
wave function of bosons can be mapped to an antisym- 
metric fermionic wave function by i^n) = 
A{xi,' - - ,xn)'^f{xi,' • •_^n) with • • • , xat) = 
Yii<j<k<N [ISl- Hence, the Hamiltonian ([1]) 
can be rewritten in terms of the fermion field operators 

Hf = ^j^dxd^¥p{x)d,i>F{x)+J^dxV{x)pF{x),{9) 

which is exactly the model describing a free fermion gas 
subjecting to the cavity periodic potential [25|. Then, 



the oscillator is formed by the particle-hole excitations 
at the edges of ±kF through the Bose- Fermi mapping, 
and the mechanical frequency is naturally related to 
the Fermi velocity: = 2kvF = 2kh7rpo/M. While 
for finite 7, we use the asymptotic expression to derive 
um = 2kvF[l - 4/7] (the dash-dotted line of Fig. [3|). 
Experimentally, when the ID gas becomes a TG gas, the 
oscillation frequency of the cavity field should follow this 
asymptotic expression and approaches cj^. Gorrespond- 
ingly, the center frequency ujs will saturate at 27r x 0.42 
MHz (Fig. EH). 

Finally, let us turn to the commensurate situation 
with Ac ~ 2/P05 where a new instability-the sine-Gordon 
transition-ma y a ppears in the strongly interacting ID 
quantum gas [19]: the superfluid ground state turns in- 
sulating in the presence of a weak commensurate periodic 
potential. It is now necessary to take account of the dis- 
crete nature of boson with m 7^ terms in Eq. JSl). This 
gives rise to a sine-Gordon type perturbation jl8| up to 
the leading term 

1 f^^ 

i^s-G = -hUoc'^cpo / dxcos[20{x) ^ Qx], (10) 

2 Jo 

where Q = 27r(po — k/ir)^ which vanishes at commen- 
surability. In the small photon numbers limit, it was 



shown that [26|, this term is renormalization irrelevant 
for K > Kc = 2 or equivalently 7 < 7c = 3.5, leaving the 
ground state a superfluid LL with the same linear exci- 
tation spectrum as the long-wavelength approximation. 
Then, one expects a well-defined oscillation of optome- 
chanical coupling. While for 7 > 7c, ^s-G becomes rel- 
evant, the system transits to an insulating Mott phase, 
where the bosonic excitations are forbidden owing to the 
energy cost, and the optomechanical oscillation vanishes 
correspondingly. Therefore, we can use the optomechan- 
ical dynamics across the critical point 7c to detect the 
sine-Gordon transition easily in experiments. 

In summary, we demonstrate that one can realize the 
effective optomechanics in the whole interacting regions 
of ID bosonic gas. This offers a new approach to to detect 
the BEC-TG gas crossover or the sine-Gordon transition 
by investigating the optomechanical coupling. These pro- 
posals are of particular significance for exploring novel 
phenomena of cavity QED and ultracold atoms. 
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